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Abstract 

We introduce a mean-field model for the potential energy landscape of a thin fluid him conhned 
between parallel substrates. The model predicts how the number of accessible basins on the 
energy landscape and, consequently, the film's ideal glass transition temperature depend on bulk 
pressure, film thickness, and the strength of the fluid-fluid and fluid-substrate interactions. The 
predictions are in qualitative agreement with the experimental trends for the kinetic glass transition 
temperature of thin films, suggesting the utility of landscape-based approaches for studying the 
behavior of confined fluids. 
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Amorphous materials confined to small dimensions play a vital role in science and technol- 
ogy. Examples include biological fluids in membranes, oil trapped in porous rocks, lubricants, 
layered composites, and thin resist films used in the fabrication of microelectronic devices. 
Many of these systems exhibit large surface-to- volume ratios, and thus their physico-chemical 
properties are influenced by boundary and finite-size effects. 1 The thermodynamic manifes- 
tations of confinement are diverse, ranging from shifts of the bulk phase boundaries to the 
creation of new phase transitions. 2 Dynamics in thin films can also differ markedly from the 
corresponding bulk materials. Most notably, the glass transition can shift to either higher 
or lower temperatures upon confinement, depending on the nature of the interactions of the 
fluid and the confining medium. 3-7 

The ability to alter the physical properties of thin films by tuning film-substrate inter- 
actions represents a tremendous opportunity for the design and fabrication of advanced 
materials. However, progress hinges on understanding processes that occur at microscopic 
and/or mesoscopic length scales. Although theoretical and computational studies continue to 
provide fundamental insights, many questions concerning the molecular origins of thin-film 
phenomena remain unresolved. 8-13 Hence, the development of a consistent and quantitative 
framework for modeling thin films is one of the outstanding challenges in the engineering 
and physical sciences. 

In this Communication, we introduce and probe the utility of an "energy-landscape" 
based approach 14 for describing the properties of thin films. Specifically, we develop a simple 
model for the topographical features of the energy landscape of a film confined between 
parallel substrates. We use this model to explore how the thermodynamic ideal glass (IG) 
transition temperature of the film (i.e., the temperature at which its configurational entropy 
vanishes) depends on bulk pressure, film thickness, and the strength of the fluid-fluid and 
fluid-substrate interactions. While the existence of an IG transition that underlies the 
laboratory glass transition of real materials remains an open and debated issue, it is a 
concept that has proven to be empirically useful and has strongly influenced modern thought 
on the glassy state. 13 We find that the predicted IG transition of our model qualitatively 
reproduces experimental trends for the kinetic glass transition of thin films. 
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FLUID FILM MODEL 



The model that we consider is the soft-sphere/mean-field (SSMF) fluid. Its potential 
energy exhibits simple scaling properties, and thus it has served as the focus of several recent 
studies (albeit, in bulk homogeneous conditions). 15,16 This fluid consists of N spherically- 
symmetric particles that interact via a soft-sphere repulsive pair potential e(ass/ r ) n hi 
addition to attractions quantified by a density-dependent mean-field form — a b p (p = N/V is 
the number density). The parameter n determines the steepness of the soft-sphere repulsion, 
and it is typically taken to be in the range 8 < n < 16. In the bulk situation, the potential 
energy per particle ip(s N , p) can be expressed in terms of the scaled coordinates of the 
particles Sj = I/ -1 / 3 ]^ (i — 1 • • • N) as 

V (s N ,p) = -a b p + y(s N K /3 , (1) 

where 7]^ = 7rcr| s p/6 is the effective packing fraction of the molecules, and y{s N ) in Eq. 
has the dimensions of energy per particle and is defined as 
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with Sij = | Si — Sj|. In the situation where two parallel substrates (each with surface area 
A) confine the particle centers to a film of volume V = AL, both the packing fraction r\ and 
the strength of the fluid-fluid attraction parameter a are functions of the film thickness L, 
and the above model is modified as 

¥ >( S N ,e L ,p) = -ap-* + y(s N )r ] n / 3 (3) 

where = <?ss/L is the dimensionless reciprocal film thickness. Here, — ^ represents 
the attraction between the fluid particles and the confining substrates, incorporated in a 
mean-field manner. Expressions for a, and rj have been derived elsewhere 17-19 and are 
given by 
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The parameter \l/o establishes the energy scale for the fluid-substrate interactions, and its 
connection to molecular parameters has been discussed previously 17,18 



POTENTIAL ENERGY LANDSCAPES OF THIN FILMS 



For given values of p and 6l, (p can be represented as a hypersurface in a 3N + 1 dimen- 
sional space — the film's potential energy landscape. Despite the multidimensional nature 
of a material's energy landscape, only a few generic features of its "rugged" topography 
have been speculated to influence the thermodynamics and dynamics of fluids. 14 " 16 Here, 
we develop a simple strategy to account for how confinement can impact these features. 
Since our primary focus is understanding the behavior of amorphous films, we restrict our 
attention to particle configurations devoid of crystalline order. 

(i) Basin enumeration function a: This function quantifies the number of distinguishable 
minima on the landscape [inherent structures). Explicitly, if the total number of (amor- 
phous) inherent structures on the landscape with well depths between <fi and <fi + d(f> is 
denoted dfl, then <r(0, e^, p) is given by the relationship dfl = C exp [Na((j>, e^, p)] d<fi, where 
C is a scale factor with dimension reciprocal energy. One of the most simple and commonly 
used approximations for the distribution of inherent structure depths in bulk fluids is the 
Gaussian function. 16 For the case of the SSMF fluid film, we propose the following simple 
generalization: 

0~oo 

where (fi m = 4> m (^L, p) and 4>oo = 4>oo{.£l, p) are the potential energy at the minimum and 
maximum values of the basin enumeration function, respectively. The above expression 
retains the functional form proposed for bulk fluids, while rendering the quantities 
and <f) m film thickness (e^) dependent. Unfortunately, the accuracy of this approximation is 
unclear at present because the landscape statistics of thin films have not been systematically 
investigated by molecular simulation. Nonetheless, we take (J5J) as a workable starting point, 
and we leave exploration of more accurate approximations for future studies. 

(ii) Mean inherent structure energy <p*: At a given (reciprocal) temperature (3 = 1/ksT, 
the fluid film will spend an overwhelming majority of its time in basins of attraction with 
inherent structures of energy (f>* = <f)*(eL,p, P), which satisfies 15 ' 16 

{j^j [0 = 0*,e L ,p]=/3. (6) 
Eq. © is exact 15,16 if the intra-basin vibrational contribution to the free energy f vi b (see 
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below) is independent of basin depth 0. For the Gaussian landscape (jSJ), we have 

<j>*{e L ,p,0) = { (7) 

/? > f3 IG 

where [3 IG = j3 IG (eL,p) locates the IG transition. 

(hi) Configurational entropy sq' This quantity is defined as sc = /cbct(0*, €l, p) and thus 
is given by 

S C (e L ,P,P) _ Woo-0m)AXoo] 2 /3</3 /G 

(iv) Ideal glass (IG) transition locus (3 IG : At temperatures below the IG transition, 
the configurational entropy of the film vanishes sc = 0, and the system is trapped in the 
amorphous basin with the lowest energy 0* = <p m : 

^%L,P) = -^- (9) 



(v) Helmholtz free energy f: In terms of the above quantities, the film's Helmholtz free 
energy possesses a simple form: 14 " 16 

f{e L , p, f3) = <f>* - s c /(3k B + f mb (10) 



IDEAL GLASS TRANSITION OF THE SSMF FLUID FILM 

Our development in the previous section briefly generalized the energy landscape formal- 
ism for bulk fluids 14 " 16 to describe fluid films. Now we quote the explicit form of the above 
functions for the SSMF fluid film. 

Using Eq. ©, (gj), and (f7|l- ([TU}l . we obtain 
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where y m and y^ represent y(s N ) of Eq. (J2J) evaluated at m and 0oo, respectively. The term 
K(/3) depends only on temperature and does not enter into our present analysis. Following 
Shell et al., ie we have modeled the vibrational contribution to the free energy of Eq. (|12j) in 
the classical harmonic approximation. 20 Other thermodynamic quantities of interest, like the 
transverse component of the pressure tensor P^e^, p, (3) = p(df /dhi p) eL ,p and the chemical 
potential p(eL, p, /3) = f + -P||P _1 of the film, follow from Eq. (fT2"j). 

To examine the confinement-induced shift of the IG transition for the SSMF fluid film, 
we consider the situation where the film is in equilibrium with a bulk fluid at pressure P b ] 
hence the "shift" we refer to is measured relative to the bulk IG transition at that pressure. 
We have determined the IG transition of the film numerically as well as by an approximate 
analytical theory. In both approaches, the IG transition of the bulk fluid is first determined 
by the condition l3l G (P b ) = (3 IG [e L = 0,pl G (P b )] (see Eq. dTTJ)), where the bulk density 
pl G (Pb) = Pb(Pb) follows from an inversion of the equation of state P b = P\\(cl = 
0,p b ,f3 b IG ). Similarly, the IG transition of the thin film (3 IG (e L ,P b ) = (3 IG [e L , p IG (e L , P h )\ is 
given by Eq. pijl . where the density of the film at its IG transition is denoted p (ei, Pb) = 
p(tL, Pb, P IG ) and is determined by the condition that the film and the bulk fluid have equal 
chemical potentials; i.e. p[eL = 0, p b (P b , P IG ), P IG ] = p, P IG ]- Numerical calculations 
require values for the parameters n, <7oo, y^, y m , a b , and \l/o- Here, we use n = 12, = 
0.5368, t/oo = 61.73e, y m = 53.22e, and a b = 16.5e<r 3 (we examine several values of the fluid- 
substrate attraction parameter ^ )- This set of parameters was chosen because it provides 
good qualitative agreement with both the liquid-state thermodynamics 16 and the predicted 
IG transition locus 21 of the bulk Lennard- Jones system. 

A perturbation approach allows us to derive the shifts in the IG transition temperature 
© = [P[ G (P b ) — f3 IG (eL, P b )]/ i3 IG (tL, Pb) to the first order in e^. To maintain brevity, we 
avoid elaborating the algebraic details that accompany these calculations, and instead quote 
the following three equivalent results: 




(13) 
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Here, Ac IG b (Pb) = —(dsc/d In f3)p b is the configurational heat capacity, K%? b (Pb) is the 
isothermal compressibility, 2T I b G (Pb)/ass = [p IG ( e L, Pb) — Pb(Pb, P ig )]/^l is the surface excess 
density, and r] b G (P b ) is the packing fraction, each evaluated at the IG transition of the bulk 
fluid, i.e. [el = 0, Pb, /3 b G (Pb)]. The thermodynamic quantities in (13) can be expressed 22 
as analytical functions of the molecular parameters of the model. Figure 1 displays the 
confinement-induced shifts in the IG transition temperature (whose physical implications 
are discussed below), demonstrating the excellent accuracy of the perturbation approach in 
capturing the full numerical results. 

The above expressions (fT3*|) shed light on the physics of confinement-induced changes in 
the IG transition temperature. Consider the first equality of Eq. (JT3J). Since Ac IG b (Pb) > 0, 
this relation reveals the following simple rule for the (small e^) shift of the IG transition. If 
confinement increases the number of basins on the landscape that the fluid can sample (hence 
increasing sc), then the IG transition temperature is depressed. Conversely, if confinement 
decreases the number of accessible basins, then the IG transition temperature is elevated. 
This result does not depend upon a specific model for the film's energy landscape. 

The second equality of Eq. ()13j) uses our model to connect the landscape-based perspective 
to physical quantities. Since d\n/3 IG /dPb < for this model, it predicts that the direction 
of the shift of the IG transition is determined by the sign and magnitude of the film's 
surface excess density 2T I b G (P b ). Large positive values of the surface excess density reduce 
the configurational entropy of the film and increase the IG transition temperature, while 
small or negative values of the surface excess have the opposite effect. While this result 
depends on our model for the energy landscape, we note that it is consistent with the 
recent theoretical predictions of McCoy and Curro. 12 The reason that these two different 
approaches arrive at similar conclusions is easy to understand. McCoy and Curro's model 
begins with the hypothesis that the confinement-induced shift in the kinetic glass transition 
is determined by how confinement affects the density of the fluid film. Similarly, as can be 
seen by Eq. fill)) , the shift in the IG transition of our mean- field model is determined by 
how confinement affects the packing fraction of the molecules in the film. 

The third equality of Eq. (|T3*j) establishes a molecular connection: strongly attracting 
walls 1 $o 'ss/ a b >> 1 elevate the IG transition while neutral or repulsive walls ^o°'ss/ a b ~ 
depress the IG transition. Fig. 1 shows the confinement-induced shifts in the IG transition 
temperature as predicted by the linear expression of Eq. (13) and the full non-linear model of 
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Eq. (11) and (12). As can be seen, the shift of the IG transition is approximately inversely 
proportional to film thickness down to molecular length scales. As has been discussed 
extensively elsewhere, these trends are in good qualitative agreement with the experimental 
shifts in the glass transition temperature of confined fluids. 3,4 ' 6 

Finally, we note that the glass transition shifts shown in Fig. 1 correspond to a fixed value 
of the bulk pressure Pb = 0. It is straightforward to employ the model and the framework 
outlined above to analyze the effects of varying the bulk pressure. Preliminary calculations 
indicate that while qualitatively similar results are seen at higher pressures, reducing the 
pressure (placing the fluid under tension) can change both the sign and the magnitude of 
shifts. 22 A comprehensive investigation of pressure effects on the glass transition is beyond 
the scope of this Communication and is deferred to a future publication. 

CONCLUSIONS 

We have examined the confinement-induced changes in the energy landscape of a simple 
statistical mechanical model of a fluid film. Our analysis here focuses on elucidating how 
the film's IG transition depends on various molecular and macroscopic parameters. The 
model predictions are qualitatively consistent with theoretical and experimental studies, 
thereby suggesting that landscape based approaches may serve to provide a framework for 
predicting the thermodynamic and dynamic properties of confined fluids. It is appropriate 
to point out that the SSMF film model, although conceptually useful, should only be viewed 
as a starting point that suggests the viability of the inherent structure formalism 14 for 
understanding the behavior of confined fluids and thin films. In fact, many films of scientific 
interest exhibit either narrow confining geometries and/or strong, directional fluid-substrate 
attractions that are not amenable to a mean-field treatment or the Gaussian landscape 
approximation. In such cases, molecular simulations can be employed to extract information 
about how confinement alters their energy landscapes, and hence their physico-chemical 
properties. We are currently pursuing research along these lines. 

We thank Scott Shell, Pablo Debenedetti, Emilia La Nave and Francesco Sciortino for 
sending us a preprint of their recent manuscript 16 , which played an important role in motivat- 
ing the present work. VG gratefully acknowledges support from National Science Foundation 
under Award Number DMR-02-04199, and the Petroleum Research Fund, administered by 
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FIGURE CAPTIONS 



FIG. 1. The confinement induced shift in the IG transition 8 as a function of for the 
SSMF fluid film. The solid lines represent the approximation of Eq. (I13|) and the symbols 
are the predictions of the full nonlinear model of Eq. and (|12|). The relative strength 
of the fluid-substrate to fluid-fluid attractions is quantified by ^o^ss/aft- The above results 
correspond to the case wherein the films are in equilibrium with the bulk SSMF fluid at 
P b = 0. 



11 




FIG. 1 - Tmskett and Ganesan, 2003 



12 



